Large current noise in nanoelectromechanical systems close to continuous mechanical instabilities 
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We investigate the current noise of nanoelectromechanical systems close to a continuous mechanical insta- 
bility. In the vicinity of the latter, the vibrational frequency of the nanomechanical system vanishes, rendering 
the system very sensitive to charge fluctuations and, hence, resulting in very large (super-Poissonian) current 
noise. Specifically, we consider a suspended single-electron transistor close to the Euler buckling instability. 
We show that such a system exhibits an exponential enhancement of the current noise when approaching the 
Euler instability which we explain in terms of telegraph noise. 

PACS numbers: 73.63.-b, 85.85.+j, 63.22.Gh 
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I. INTRODUCTION 

Nanoelectromechanical systems (NEMS) in which a 
nanomechanical resonator is coupled to electronic degrees of 
freedom 13 through, e.g., a single-electron transistor (SET), 
show spectacular effects stemming from the coupling of the 
mechanical part of the device to the electronic charge. These 
effects arise due to the reduced size of the nanoresonator, 
so that the backaction of the mechanical degree of freedom 
on the SET can have significant consequences for the trans- 
port properties. A prominent example is the low-bias current 
blockade that occurs in the Coulomb blockade regime when 
the nanoresonator is capacitively coupled to the SET. 4 ' 5 The 
presence of an extra electron with charge — e < on the cen- 
tral island forming a quantum dot on the suspended vibrating 
structure induces an electrostatic force F c on the resonator, 
shifting the equilibrium position of the latter by an amount 
F c /k, with fc the spring constant of the oscillator (see Fig. 1). 
This induces a shift of the gate voltage V g ~ F^/ek seen by 
the SET, and, hence, a blockade of the current through the 
device for bias voltages V < F^/ek. This phenomenon is 
the classical counterpart of the Franck-Condon blockade in 
molecular devices 6 7 that has recently been observed in car- 
bon nanotube-based resonators for high-energy longitudinal 
stretching modes. 8 For classical nanoresonators, the current 
blockade has, to the best of our knowledge, never been ob- 
served experimentally due to the relatively weak electrome- 
chanical coupling F c to the low-energy bending modes of the 
suspended structure, although a precursor of this effect has 
been lately reported in the literature. 9 ' 10 

It has been recently shown 1112 how one can enhance the 
classical current blockade by orders of magnitude by exploit- 
ing the well-known Euler buckling instability. 13 Indeed, the 
spring constant k (or equivalently, the vibrational frequency 
of the fundamental bending mode w) tends to zero when one 
brings the nanoresonator to the buckling instability with the 
help of a lateral compression force F (see Fig. 1). Thus, 
the energy scale F^/k at which the current blockade occurs 
dramatically increases, rendering this phenomenon potentially 



FIG. 1. (Color online) Sketch of a suspended doubly clamped 
nanobeam forming a quantum dot (solid black line) connected via 
tunnel barriers to source and drain electrodes held at chemical po- 
tentials /il and fin by the bias voltage V. The lateral force F com- 
presses the nanobeam and induces the buckling instability. The beam 
is capacitively coupled to a metallic electrode kept at a gate voltage 
V g . This induces a stochastic force F c that attracts the beam towards 
the gate electrode whenever the quantum dot is charged (dashed red 
line), inducing fluctuations of the nanobeam's deflection and, in turn, 
fluctuations of the current through the device. 



observable in future experiments. 

It is the purpose of the present paper to investigate the cur- 
rent noise in the vicinity of a mechanical instability, using 
the Euler buckling instability as a paradigmatic model. We 
find that the current noise (which contains valuable informa- 
tion about the dynamics of the nanomechanical system 14-20 ) 
is strongly enhanced in the vicinity of the Euler instabil- 
ity. The underlying source of noise arises from the stochas- 
tic nature of the charge transfer processes. These are pro- 
ducing a current-induced stochastic force acting on the me- 
chanical degrees of freedom, consisting of current-induced 
(conservative and nonconservative) averages as well as a 
fluctuating force. 4 ' 5 ' 1 X ' 12 ' 16 » 21 - 26 Hence, the deflection of the 
nanotube exhibits a Langevin dynamics which, due to the 
backaction of the nanoresonator on the SET, produces large 
super-Poissonian current noise. 5 ' 14 ' 16 ~ 20 This effect is particu- 
larly strong close to the Euler buckling instability, where the 
nanoresonator becomes extremely soft (fc — > 0). 
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Our theoretical approach employs a nonequilibrium Born- 
Oppenheimer approximation, 4 ' 516 ' 2125 which exploits the 
separation of timescales between fast electronic and slow 
mechanical degrees of freedom. This leads to an effec- 
tive stochastic description of the nanoresonator in terms of 
a Langevin equation. This approach becomes asymptoti- 
cally exact in the vicinity of the mechanical instability where 
u -► O. 11 ' 12 - 26 

The paper is organized as follows: Our model and the ef- 
fective Langevin description of the nanoresonator deflection 
is presented in Sec. II. In Sec. Ill, we briefly recall the main 
results of Ref. 1 1 for the current-voltage characteristics of the 
system that are essential for the understanding of our numer- 
ical investigation of the current noise presented in Sec. IV. In 
Sec. V, we detail the role played by thermal fluctuations on the 
current noise and propose an analytical model based on tele- 
graph noise that reproduces most of our numerical findings. 
We present in Sec. VI the role played by the full nonequilib- 
rium fluctuations on the noise before we conclude in Sec. VII. 



H. MODEL 

The model we adopt is the same as in Ref. 1 1 and we briefly 
recall it here for the convenience of the reader. The setup (see 
Fig. 1) consists of a quantum dot embedded in a suspended 
nanobeam connected via tunnel barriers to source and drain 
leads. A lateral compression force F exerted on the nanobeam 
brings it to the Euler buckling instability when F approaches 
the critical force F c . The gate electrode induces an electrome- 
chanical coupling between the bending modes of the tube and 
the charge state of the dot proportional to the electrostatic 
force F c . The Hamiltonian of the system reads 

H = -ffvib + HsET + He, (1) 

where i/ v ib describes the oscillating modes of the nanobeam, 
.ffsET the single-electron transistor, and H c the coupling be- 
tween vibrational modes and electronic degrees of freedom. 

At the Euler buckling instability (F = F c ), the frequency 
of the fundamental bending mode 

Tp 

ui = UqV^S, 5 = — -1 (2) 

vanishes, 13 ujq being the frequency of that mode for F = 0, 
while all higher energy modes have a finite frequency. At suf- 
ficiently low temperatures, we thus only consider the funda- 
mental bending mode and write 11 ' 12 ' 26 " 30 

for the vibrational part of the Hamiltonian (1). Here, X is 
the deflection of the tube and P its associated canonical mo- 
mentum with effective mass m. In Eq. (3), the quartic term 
proportional to a > ensures the stability of the system for 
F > F c (cj 2 < 0), where the beam buckles into one of the 
two metastable positions at X = ±^J —muj 2 / a. For F < F c 
(w 2 > 0), the beam remains flat. 



We model the SET by a single-level quantum dot with or- 
bital energy ed connected via tunnel barriers to left (L) and 
right (R) leads held at chemical potential /il and /xr, respec- 
tively. The SET Hamiltonian reads -ffsET = -Hdot + -fficads + 
-Htun, where the dot Hamiltonian Hdot = ( e d — eV g )n d + 
Unwind — l)/2 is expressed in terms of n d = d^d, d an- 
nihilating an electron on the dot. Here, V s is the (effective) 
applied gate voltage. The intradot Coulomb repulsion is de- 
noted by U and is assumed to be the largest energy scale in 
the problem, thus preventing double occupancy of the dot. 
The Hamiltonian for the (spinless) electrons with energy e k 
and momentum k in the two leads (annihilated by the operator 
c ka , a = L, R) is written as H lcads = Y,ka( € k ~ Va)c\ a c ka . 
Finally, tunneling between dot and leads is accounted for by 
the Hamiltonian H tun = ^ ka {t a <J ka d+\\.c), with t a the tun- 
neling amplitude between the dot and lead a. In the remain- 
der of this paper, we assume the temperature T to be much 
larger than the tunneling-induced width T = J2 a of the 
dot orbital (sequential tunneling regime). This weak-coupling 
assumption should not qualitatively change our results for the 
low-frequency current noise, as it is the case for the I-V char- 
acteristics which is qualitatively the same in the sequential 11 
and cotunneling 12 regimes. Moreover, we consider for sim- 
plicity symmetric voltage drops (/zl, = — /xr = eV/2) and 
symmetric coupling to the leads (IY = Tr = T/2). 

The coupling Hamiltonian between vibrational and elec- 
tronic degrees of freedom entering Eq. (1), 

H c = F c Xn d , (4) 

arises due to the capacitive coupling between the gate elec- 
trode and the nanobeam when the latter is charged with one 
extra electron. 9 ' 10 ' 31 Since the dot occupation n d is a stochas- 
tic variable (taking values and 1 in our model), the coupling 
(4) produces a random electrostatic force with magnitude F c 
on the nanobeam (see Fig. 1). 

As uj — » close to the buckling instability, the 
oscillator becomes classical (hut -C fceT) and slow 
compared to the electronic degrees of freedom (u <C 
T). This justifies a nonequilibrium Born-Oppenheimer 
approximation, 11 ' 12 ' 16 ' 21 " 26 in which the vibrational dynamics 
is characterized by a Langevin process with white noise, 

^ + [7(z)+7e]^=/rf(&)+£foT). (5) 

Here and in what follows, we use reduced units x — X/£, 
p = P/mu>o£ and r = ujQt in terms of the polaron shift £ = 
F e /mujQ and the relevant energy scale of the problem = 
F c £ (for more details, see Ref. 11). In Eq. (5), the effective 
force acting on the nanobeam, 

f e s(x) = 8x - ax 3 - n (x), (6) 

with a — a£ 4 /E^, arises (i) from the bare vibrational Hamil- 
tonian (3) and (ii) from the coupling between vibrational and 
electronic degrees of freedom (4). This current-induced force 
is proportional to the occupation of the dot for fixed x, which, 
in the sequential tunneling regime (HT -C k^T), is given by 

no(x) = ^[h(x)+f K (x)}, (7) 
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with 



(8) 



the Fermi function in the left and right leads, respectively. 
Here we introduced a reduced bias voltage v = eV/E^, gate 
voltage v g — (eV g — e^)/E^, and temperature T = k^T/E^. 
The charge fluctuations on the quantum dot induce a fluctuat- 
ing force £(x,t) in the Langevin equation (5), with average 
(£(ie,t)) = and white-noise correlator t)£(x, t')) = 
\d(x) + 2j c T]5(t — t'). Here, the current-induced fluctuation 
is given by 4,11 



d{x) 



2u! 



n (x) [1 - n (x)} 



(9) 



and the extrinsic damping constant 7 C accounts for the finite 
quality factor Q = 1 /j e of the nanoresonator. Finally, retar- 
dation effects in the response of the resonator to the current 
flow lead to a current-induced dissipative force with friction 
coefficient 411 



1 ox 



(10) 



The Langevin equation (5) is equivalent to the Fokker- 
Planck equation 32 



d_ 



V{x,p,r) = CP{x,p,r) 



(11) 



for the probability distribution V (x, p, r) that the nanobeam is 
at phase-space point (x,p) at time r. In Eq. (11), the Fokker- 
Planck operator is defined as 



C 



P dx- fcs{x) l P + [lix) 



dx 
d{x) 



1 d 
7c] dp P 



dp 2 



(12) 



Solving the Fokker-Planck equation (11) [or equivalently the 
Langevin equation (5)] gives access to both the dynamics of 
the vibrational mode of the nanoresonator and the resulting 
transport properties of the device, such as its I-V characteris- 
tics (see Sec. Ill) and its noise power spectrum (see Sees. IV, 
V and VI). 



for fixed deflection x over the stationary solution V s t of the 
Fokker-Planck equation (11). 

The classical current blockade phenomenon 4 ' 11 can be un- 
derstood in terms of the effective potential 



. . 5x 2 ax 4 
v eS {x) = — — + — 



|ln(/ L (z)/ R (z)) 



(15) 



associated to the effective force (6). The effective potential is 
shown in Fig. 2 for T = together with the average occupa- 
tion of the dot (7) and the quasistationary current (14) for com- 
pression forces corresponding to the beam far below [Figs. 
2(a)-2(c)] and at the Euler instability [Figs. 2(d)-2(f)]. 33 In 
both cases, the most stable minima of v e g(x) correspond, for 
bias voltages v smaller than the gap 11 



A, 



-1/25, 
1/46, 

2 1/3 -l / 3 

5V3 



& - 

AV3, 



<5»a 1 /3 j 
\5\ « a 1 / 3 , 



(16) 



to an average occupation n Q (x) = or 1, i.e., the system 
is not conducting ["blocked" minima for which I(x) = 0, 
cf. Figs. 2(b) and 2(e)]. For v > A„, the most stable mini- 
mum corresponds to no(x) = 1/2 and the current can flow 
["conducting" minimum corresponding to I(x) ^ 0], At the 
threshold v = A v , the three minima are metastable, leading to 
a current-induced instability of the system. Since for relevant 
experimental parameters, a <C l, 9_n the gap (16) is maximal 
at the instability where 6 = [F = F c , cf. Eq. (2)] and is or- 
ders of magnitude larger than for 5 — — 1 (F — 0). The gaps 
of Eq. (16) are obtained for gate voltages v g = i 1 ™ 1 ", with 11 



1/25, -5>a 1 / 3 , 
-1/48- y/5j&, 5>a 1 / 3 , 



(17) 



At finite temperatures, the effective potential (15) qualita- 
tively changes (see Fig. 3): the barriers separating the min- 
ima are lowered as one increases the temperature, leaving the 
effective potential with a single minimum for large enough 
temperatures. As detailed in Ref. 11, the current blockade 
becomes less pronounced as the temperature increases [see 
Fig. 8(a) below]. Moreover, taking into account the current- 
induced fluctuations (9) and dissipation (10) in the Langevin 
dynamics (5) further reduces the current blockade [see Fig. 
11(a) below]. 



III. CURRENT BLOCKADE 

Due to the separation of timescales between slow vibra- 
tional motion and fast electronic dynamics, the average cur- 
rent 



1 = J dxdpV st {x,p)l{x) (13) 



is obtained by averaging the quasistationary current 



l{x) 



eT 



lh(x)-f R (x)} 



(14) 



IV. CURRENT NOISE AT THE EULER BUCKLING 
INSTABILITY 

We start by discussing the two contributions to the current 
noise, i.e., the usual shot noise and the mechanically-induced 
noise. It is the latter contribution to the noise which we find to 
be dramatically enhanced close to the mechanical instability. 

The noise power spectrum is defined as 34 



s{n) = 2 



+oo 



dt e lUt (AI(t +t)AI(t )), 



(18) 
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FIG. 2. (Color online) Zero temperature (a,d) average occupation of the dot for fixed x [Eq. (7)], (b,e) quasistationary current [Eq. (14)] and 
(c,f) effective potential [Eq. (15)], (a,b,c) far below the Euler instability (—5 ^> a 1//3 ) and (d,e,f) at the instability (S = 0) for increasing bias 
voltages and for a gate voltage v & = v™ ln [cf. Eq. (17)]. 
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FIG. 3. (Color online) For increasing temperature, (a,d) average occupation of the dot for fixed x, (b,e) quasistationary current and (c,f) 
effective potential, (a,b,c) far below the Euler instability (6 = — 1) and (d,e,f) at the instability (5 — 0). In the figure, v — A v , u g = f™ 111 and 
a = 10~ 6 . 



where AI(t) = I(t) — (I) denotes the time-dependent fluctu- 
ations of the current operator. In Eq. (18), the brackets (. . .) 
indicate an ensemble average or, equivalently, an average over 
the initial time to. Due to the separation of timescales be- 
tween fast electronic dynamics and slow vibrational motion 
(lu <C r), one can identify two additive contributions to the 
noise power spectrum (18), S = S sh + S^. 5 - 16 - 20 - 21 The first 
one, S^, corresponds to the (thermal) Nyquist-Johnson and 
shot noise and is discussed in Appendix A. The second con- 



tribution to Eq. (18), the mechanically-induced noise S m (re- 
ferred to as "mechanical noise" in the sequel), is induced by 
the fluctuations of the nanobeam deflection x[t). These occur 
on a much longer timescale (of the order of 1/w) than the shot 
noise (the corresponding current-current correlator decaying 
in that case on the short timescale 1/r). The mechanically- 
induced noise therefore dominates the noise power spectrum 
at low frequencies, and can exceed the shot noise by orders of 
magnitude . 5 ' 1 6 ' 20 ' 2 1 
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The mechanical noise reads 

S m (fi) = 2 J -~ e inr ^° J dxdpdx dp AI(a;) 

x "P(x,p, T\xo,po,T )Al(xo)V s t(xo,po) (19) 



where V(x,p, t\xq,Pq, Tq) is the conditional probability that 
the nanobeam is at phase-space point (x,p) at time r, pro- 
vided it was at (xo,po) at time To = 0. In Eq. (19), 
Al(x) = I(x) — I is the quasistationary current fluctuation, 
withZ(a;) and / given by Eqs. (14) and (13), respectively. The 
conditional probability "P(:r,p, t\xq,pq, to) can be obtained 
from the time-dependent solution of the Fokker-Planck equa- 
tion (11) with the initial condition V(x,p : To\xo,po, r o) = 
S(x — Xo)S(p — po). Equation (19) can then be re-expressed 
by exploiting the above initial condition and performing the 
Laplace transform of the Fokker-Planck equation (11). This 
procedure yields 5 



Sm(«) 



4 

w 



dxdp AI(x) 



C[Al(x)V st (x,p)], (20) 



with C the Fokker-Planck operator defined in Eq. (12). 35 

In the form of Eq. (20), the mechanical noise can be 
straightforwardly computed numerically, since it only re- 
quires the knowledge of the stationary probability distribution 
■p s t corresponding to the Fokker-Planck equation (11), as it is 
the case for the average current (13) (for details, see Ref. 5). 36 
Alternatively, one can solve for the time-dependent solution 
x(t) of the Langevin equation (5) using standard techniques 
for stochastic differential equations. 37 The average current and 
noise are then computed by performing the time averages of 
the quasistationary current (14) and the current-current cor- 
relator entering Eq. (18), respectively. We have checked nu- 
merically that this approach yields the same results as the ones 
based on the stationary solution of the Fokker-Planck equation 
(11) [see Eqs. (13) and (20)]. However, although more phys- 
ically transparent, the method based on Eq. (5) requires long 
simulation times as well as sampling. We thus use the other 
method for all the numerical results presented in the sequel of 
the paper [except in Fig. 5, where we explicitly simulate the 
time-dependent deflection of the nanobeam from the Langevin 
equation (5)]. 

In Fig. 4, we present our numerical results for the Fano 
factor F m — S m (0)/2e\I\, where the zero-frequency noise 
S m (0) and the average current /, Eqs. (20) and (13), respec- 
tively, are computed for typical parameters as a function of 
the reduced force 5. In Fig. 4, the bias and gate voltages cor- 
respond to the apex of the Coulomb diamond [v = A v and 
V g = f™ m , cf. Eqs. (16) and (17), respectively]. As envi- 
sioned above, there is a dramatic increase of the mechanically- 
induced current noise in the vicinity of the Euler buckling in- 
stability (8 w in Fig. 4) as compared to the noise far away 
from the instability. Moreover, the Fano factor can take, de- 
pending on the compression force 6, super-Poissonian values 
(F m > 1) that are well above the shot noise contribution, 
F s h = 1/2 (see Appendix A). 
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FIG. 4. (Color online) Fano factor as a function of the reduced com- 



pression force S [cf. Eq. (2)]. In the figure, 7 C = ujq/T 
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T = 3, a 
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A„ and v g 



The red circles corre- 



spond to the data points, while the dotted line serves as a guide to the 
eye. 



The results of Fig. 4 can essentially be understood in terms 
of telegraph noise in the effective potential (15) (see also Figs. 
2 and 3). 38 Indeed, unlike the energy gap (16) which varies al- 
gebraically with the force S as ~ 1 / j <5| , the numerical results 
of Fig. 4 indicate that the noise (or Fano factor) depends expo- 
nentially on l/\8\ (notice the logarithmic scale in Fig. 4), sug- 
gesting telegraph noise. As the compression force increases 
towards the instability at 5 = 0, the height of the barriers sep- 
arating the three metastable minima (for v = A v ) grows as 
the gap (16), such that the waiting time of the system in one 
of these minima increases exponentially. Thus, the probabil- 
ity for the system to switch to another minimum is drastically 
reduced, subsequently increasing the telegraph noise. As the 
height of the potential barriers near the Euler instability is very 
large, scaling as 1/<5 1/3 with a < 1 [see Fig. 2(f)], and the 
energy gap (16) is maximal at the instability, this leads to a 
current noise which is also maximal at the Euler instability. 

This interpretation is confirmed by Fig. 5, which shows the 
result of a simulation 37 of the deflection of the nanobeam x 
as a function of time (see blue lines in the figure) as obtained 
from the Langevin equation (5) for the same parameters as in 
Fig. 4. We also show the resulting quasistationary current (14) 
as a function of time by red lines in Fig. 5. Far from the in- 
stability [Figs. 5(a) and 5(e)], the dynamics of the nanobeam 
follows qualitatively the behavior of a Brownian particle in a 
harmonic potential. Indeed, for the temperature used in Figs. 
4 and 5, the effective potential (15) far from the instability has 
a single minimum [see also Fig. 3(c)]. For temperatures T 
which are large compared to the gap (16), the current shown 
in Figs. 5(a) and 5(e) switches rapidly between values which 
are small as compared to the maximal current eT/4 [cf. Fig. 
3(b)]. Hence, the resulting Fano factor is relatively small (cf. 
Fig. 4). As one approaches the Euler instability from below 
[Fig. 5(b)] or above [Fig. 5(d)], the dynamics of the nanobeam 
becomes slower. Then the behavior of the current as a func- 
tion of time starts to resemble telegraph noise, as the effective 
potential starts developing metastable minima for this value 
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FIG. 5. (Color online) Deflection x (in blue) and quasistationary 
current T (in red) as a function of time simulated by the Langevin 
equation (5). In the figure, the compression force increases from (a) 
to (e): (a) 5 = -1, (b) S = -0.05, (c) 5 = 0, (d) 5 = 0.05, and (e) 
(5 = 1. The parameters used in the figure are the same as in Fig. 4. 



of the temperature as compared to the energy gap (16). At 
the instability [Fig. 5(c)], the dynamics of the nanobeam be- 
comes very slow, and the behavior of the current as a function 
of time is completely stochastic and uncorrelated, with long 
waiting times between vanishing and maximal current. The 
corresponding Fano factor is thus extremely large and super- 
Poissonian (F m > 10 3 in Fig. 4), and much larger than far 
from the Euler instability. 

In order to understand the features of our main numerical 
results presented in Figs. 4 and 5 in more detail, we will first 



consider the role played by thermal fluctuations alone (Sec. 
V), while the full nonequilibrium dynamics is investigated fur- 
ther in Sec. VI. 



V. THERMAL FLUCTUATIONS 

In this section, we consider the fully adiabatic limit 
wo /r — » 0. We neglect the current-induced fluctuations and 
dissipation in the Fokker-Planck equation (11) [cf. Eqs. (9) 
and (10)] and consider the role played by thermal fluctuations 
alone. We start in Sec. V A with a simplified analytical model 
based on standard telegraph noise. In Sec. V B we substantiate 
our findings by evaluating Eq. (20) numerically. 



A. Telegraph noise 

We present here an analytical estimate of the noise power 
spectrum. Our simplified model relies on thermally-induced 
telegraph noise 39 and on an estimate of the escape rates based 
on Kramers reaction rate theory. 32 ' 40 ' 41 

In what follows, we work in the low-temperature regime 
T <C A u , with the gap A t , given in Eq. (16). Moreover, 
we focus on the case where the nanobeam is far below the 
Euler instability (—5 a 1 / 3 ), as the results presented be- 
low should stay qualitatively the same for larger compression 
forces. We thus approximate the effective potential (15) by its 
zero-temperature counterpart, which is shown in Fig. 2(c) for 
a gate voltage v g = v™ m [cf. Eq. (17)]. As one can see from 
Fig. 2(c), the effective potential has three metastable minima 
for < v < 2A V : two of them are equivalent (located sym- 
metrically at x\ = l/J and xq = about the line x — 1/2(5) 
and correspond to a state in which the current vanishes [see 
Fig. 2(b)], while the one at 2^/2 = 1/2(5 corresponds to a 
current-carrying state. This suggests to write a rate equa- 
tion for the probabilities P c and P^ that the system is in a 
conducting or blocked state, respectively. Denoting Tin = 

(-Tx — >X l/2 anC l Tout = Y Xl/2 ^ Xl + I" 

the transition rates in and out of the conducting state (T Xi - ¥Xj 
is here the transition rate from the minimum located at x; to 



the one at xa), we have P c 



-Ph 



-T n „tP r . + FinR, 



Following Ref. 39, the average current and the noise power 
spectrum are readily obtained from the above rate equation. 
They read 



I = 



e.T 



4 r 



and 



Sm(fi) 



■Tin To 



4 r in + r out ^ + (r. 



(21) 



(22) 



respectively. Notice that for bias voltages v ^ 2 A.„, the effec- 
tive potential (15) has a single minimum [see Fig. 2(c)], and 
the telegraph noise model presented above does not apply. In- 
stead, the system's dynamics is characterized in that case by 
standard Brownian noise. 
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As detailed in Appendix B, the transition rates entering 
Eqs. (21) and (22) can be easily calculated using Kramers 
theory. 3240 ' 41 Incorporating Eq. (B2) in Eq. (21), we find for 
the average current the approximate expression 



I = 



eT 

T 



1 + 2 cxp 



4T 



1 



v 

a7 



(23) 



valid for v < 2A V . Using Eqs. (22) and (B2), we find for the 
zero-frequency noise 42 



S m (0) = 




(24) 



where the function h(z) is defined in Eq. (B3). 

At a bias voltage corresponding to the energy gap (16) (v = 
A v ), the results of Eqs. (23) and (24) simplify to yield the 
Fano factor 



ier 

3w 



7e 



T 

-r— exp i — =- , . 
A„ V16T/ 



(25) 



Although this result is based on a simplified model and despite 
the fact that it does not include the full nonequilibrium dynam- 
ics of the nanoresonator induced by the charge fluctuations on 
the dot, it qualitatively captures our main finding depicted in 
Fig. 4. Indeed, as one approaches the Euler instability from 
below, the gap A^ increases algebraically as ~ 1/|<5| [cf. Eq. 
(16)], resulting in an exponential increase of the Fano factor. 

The results of Eqs. (23), (24) and (25) also apply for com- 
pression forces far above the Euler instability (i5 a 1 / 3 ). 33 
Since the energy gap A„ is, in that case, half of the gap far 
below the instability [—5 a 1 / 3 , cf. Eq. (16)], this explains 
the asymmetry of the Fano factor for negative and positive S in 
Fig. 4. In the vicinity of the buckling instability (|<5| <C a 1 / 3 ), 
the exponential dependence of the Fano factor as a function 
of the gap (16) (which here scales with a < 1 as 1 /a 1 / 3 ) 
should stay qualitatively the same. This results in a Fano fac- 
tor which saturates at its maximal value at the Euler instability 
(cf. Fig. 4). 

The results of Eqs. (23) and (24) for the average current 
and the zero-frequency noise are shown in Figs. 6(a) and 6(b), 
respectively. As one can see from Fig. 6(b), our analytical re- 
sults capture the following trends for the mechanical noise: (i) 
it only depends on the compression force 5 through the ratio 
T / A v , (ii) the noise is maximal close to v — A v (a feature 
which has also been found in Ref. 5 in the case of fully coher- 
ent transport) and its maximum shifts towards higher bias volt- 
ages when one increases the temperature, and (iii) the noise is 
inversely proportional to the extrinsic damping constant (i.e., 
proportional to the quality factor). In Sec. V B, these features 
based on our simplified telegraph-noise model will be con- 
firmed and discussed further in the context of numerical cal- 
culations based on Eq. (20). 
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FIG. 6. (Color online) (a) Current (23) and (b) zero-frequency noise 
(24) as a function of bias voltage for increasing temperature as ob- 
tained from the telegraph-noise model. 



B. Numerical results 

In Fig. 7, we present numerical results for the average cur- 
rent [Fig. 7(a)] and the zero-frequency noise [Fig. 7(b)] as a 
function of bias voltage, for a gate voltage corresponding to 
the apex of the Coulomb diamond, see Eq. (17). As one can 



to 




FIG. 7. (Color online) (a) Current and (b) zero-frequency noise as 
a function of bias voltage for increasing values of the reduced com- 
pression force 8. In the figure, v e — w™ m , T — A„/10, 7 C = 10~ 2 , 
cuo/r — (no current-induced fluctuations), and a — 10~ 6 for all 
lines except for the thick blue dashed one, where a = 10 . 
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FIG. 9. (Color online) Zero-frequency noise as a function of bias 
voltage for increasing values of the extrinsic damping constant (or 



inverse quality factor) 7 e . In the figure, S 



-1, v K 



T 



A„/10, tJo/r = 0, and a = 10 6 . Inset: Same as the main figure 
for 6 = 0. 



FIG. 8. (Color online) (a) Current and (b) zero-frequency noise as a 
function of bias voltage for increasing temperature T = 0.03 (solid 
line), 0.04 (dashed line), 0.06 (dotted line), 0.08 (dash-dotted line). 
In the figure, 5 = -1, u g = v™ in , 7 C = 10" 2 , iv /F = 0, and 
a = 10~ 6 . Insets: Same as the main figure for 6 = 0. In the insets, 
f = 2 (solid line), 2.5 (dashed line), 3 (dotted line), 3.5 (dash-dotted 
line). 



see from the figure, the noise has qualitatively the same behav- 
ior far from the Euler instability [\5\ a 1 / 3 , see thin black 
lines in Fig. 7(b)] and at the instability [6 — 0, thick lines in 
Fig. 7(b)]. Remarkably, the noise [as well as the current, see 
Fig. 7(a) and Ref. 11], once plotted as a function of v/A v , 
and for the same value of T/A v , is (almost) quantitatively the 
same for any S far from the instability [see thin black lines in 
Fig. 7(b)]. The overall behavior of the current and noise in 
Fig. 7 when the system is far from the Euler instability (see 
thin solid black lines in Fig. 7) can be understood in terms 
of the telegraph-noise model detailed previously. Indeed, Eq. 
(24) shows that for T <C A„, the noise is exponentially sensi- 
tive to the ratio A v /T only, thus explaining the scaling in Fig. 
7(b). In contrast, the above scaling does not apply at the me- 
chanical instability [see thick red and dashed blue lines in Fig. 
7(b)]. There, the noise is larger for smaller values of the an- 
harmonicity parameter a, the latter determining the strength 
of the quartic correction in the effective potential (15). 

In Fig. 8, we present the temperature dependence of both 
the current and the zero-frequency noise far below the Euler 
instability [Figs. 8(a) and 8(b), respectively] and at the insta- 
bility (see insets in Fig. 8). As one can see from the figure, the 
behavior of these two quantities is qualitatively the same far 
from and at the buckling instability. As temperature increases, 
the current blockade becomes less pronounced [see Fig. 8(a)] 
as the system can explore more phase-space due to thermal 
fluctuations (for more details, see Ref. 11). Moreover, the 
maximum of the zero-frequency noise decreases with temper- 
ature and gets shifted to larger values of the bias voltage [see 



Fig. 8(b)], a phenomenon which is captured by our analytical 
estimate of the noise in Sec. VA (see Fig. 6). As tempera- 
ture increases, the probability to jump out of one minimum 
of the effective potential (cf. Fig. 2) increases exponentially, 
such that the associated telegraph noise decreases. 

In Fig. 9, we show numerical results for the zero-frequency 
noise for various values of the extrinsic damping constant 
7c = X/Q {ox inverse quality factor) far below the instabil- 
ity (Fig. 9) and at the Euler instability (inset in Fig. 9). As 
one can see from the main figure, the mechanically-induced 
noise scales almost perfectly as 5 m (0) ~ 7" 1 for compres- 
sion forces far below the instability. This is characteristic 
of telegraph noise in the weak friction limit, 40 ' 41 where the 
system's energy varies only slowly with time [cf. Appendix 
B and Eq. (24)]. Interestingly, in the case of molecular de- 
vices, the noise is also much larger for unequilibrated (high- 
er) vibrons. 6 ' 7 At the instability (see the inset in Fig. 9), the 
scaling with 7 e is only approximate. 43 

We conclude this section by computing the frequency de- 
pendence of the mechanical noise (20) far below (Fig. 10) and 
at the Euler instability (inset in Fig. 10). Far from the insta- 
bility, the frequency dependence of the noise shows a l// 2 
dependence, typical of telegraph noise [cf. Eq. (22)]. Notice 
that the width of the Lorentzian shape of 5 m (f2) depends on 
the bias voltage through the transition rates for the system to 
enter or leave the conducting minimum corresponding to an 
average occupation of the dot of 1/2, see Eqs. (22) and (B2). 
At the Euler buckling instability, 5 m (f2) also follows a l// 2 
behavior for £1 <C Wq, although additional structures appear 
for larger frequencies and for certain bias voltages. 44 



VI. NONEQUILIBRIUM FLUCTUATIONS 

We now investigate the mechanical noise in the presence of 
the current-induced fluctuations and dissipation, Eqs. (9) and 
(10). 
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FIG. 10. (Color online) Frequency dependence of the mechanical 
noise for various values of the bias voltage. In the figure, 8 = — 1, 

v g = vf n , f = A„/10, a = 10~ 6 , 7 C = 10~ 2 , and u) /T = 0. 
(Inset) Same as the main figure for 5 = 0. 



Our numerical results for the current and the zero- 
frequency noise are shown in Figs. 11(a) and 11(b), respec- 
tively, when the nanobeam is at the Euler instability (the in- 
sets in Fig. 11 consider the case 5 = —1). As one can see 
from Fig. 11, the effect of an increasing adiabaticity parame- 
ter ujq /r, which controls the strength of the current-induced 
fluctuations (9) and dissipation (10), is qualitatively similar 
to the effect of an increasing temperature (cf. Fig. 8). In- 
deed, the overall noise level is reduced and the maximum of 
the noise is shifted towards larger bias voltages for increasing 
cjo/r [see Fig. 11(b)]. Moreover, the current blockade gets 
less pronounced for increasing ujo/T [see Fig. 11(a)]. 

The results of Fig. 1 1 can be qualitatively understood by 
defining an effective temperature 11 ' 12 

- (d)/2 + 7c f 

Jeff = , . (26) 

(7) + 7c 

in analogy with the fluctuation-dissipation theorem. 32 Here, 
(d) and (7) are the averages over phase-space of the current- 
induced fluctuations and dissipation, Eqs. (9) and (10), respec- 
tively. As shown in Ref. 11, the effective temperature (26) is 
approximately given by T c ff — T+ "^ r ®(v— A v ), with 9(z) 
the Heaviside step function. This explains why both the cur- 
rent and the zero-frequency noise are quite insensitive to the 
ratio wo /T for v < A v and are similar to the case ujq /Y = 0, 
i.e., the fully adiabatic limit. On the contrary, for v > A„, 
the effective temperature increases with increasing ojq /T, ex- 
plaining the similarity of the behavior of the current and noise 
in Figs. 1 1 and 8. 

We conclude this section by noticing that we have numer- 
ically checked that the frequency dependence of the noise 
power spectrum also follows a l// 2 behavior when one takes 
into account current-induced fluctuations, similar to Fig. 10. 
This confirms that the noise is dominated by a telegraph noise 
at low-enough temperatures even in presence of nonequilib- 
rium fluctuations. 




12 3 4 



v/A v 

FIG. 11. (Color online) (a) Current and (b) zero-frequency noise as 
a function of bias voltage for increasing values of the adiabaticity 
parameter uio/T. In the figure, 6 = 0, v e = i>™ m , 7c = 10 -2 , 
T/A v =0.1 and a = 10~ 6 . (Inset) Same as the main figure with 

S=-l. 



VII. CONCLUSION 

We have investigated the current noise in nanoelectrome- 
chanical systems close to a continuous mechanical instability. 
Specifically, we have considered the paradigmatic Euler buck- 
ling instability in suspended single-electron transistors which 
are capacitively coupled to a gate electrode. We have pre- 
dicted a drastic enhancement of the current noise when the 
nanobeam supporting the quantum dot is brought to the Euler 
instability, resulting in very large Fano factors that are well 
above the Poisson limit. This exponential enhancement at the 
buckling instability is directly related to the (algebraic) en- 
hancement of the current blockade predicted in Ref. 1 1 . We 
developed a rather detailed picture of the underlying physics 
in terms of a telegraph-noise model. While such large Fano 
factors may make observation of the low-bias current block- 
ade more challenging, the large noise levels predicted in this 
work would serve also as a clear experimental signature of the 
interplay between electronic and mechanical degrees of free- 
dom in NEMS (e.g., carbon nanotubes 45 ) close to continuous 
mechanical instabilities. 
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FIG. 12. (Color online) Zero-frequency shot noise as a function of 
bias v and gate voltage v g (a) for vanishing compression force F = 
and (b) at the Euler buckling instability where F — F c . The param- 
eters used in the figure are T = 10, a = 10~ 6 , wo/r = 10~ 2 , and 
7c = l/Q = 10 -2 . Color scale: white and red (light gray) regions 
correspond to S s h(0) = and e 2 T/4, respectively. 



takes the value 5 s h(0) = e 2 T/4 = e\J\ in the conducting 
regions of the v-v g plane, with a corresponding Fano factor 
F s h = 5 s h(0)/2e|/| of 1/2, typical for a single-level quan- 
tum dot symmetrically coupled to source and drain leads in 
the sequential-tunneling regime. 34 The behavior of the mean- 
field zero-frequency shot noise as a function of the reduced 
compression force 6 can thus readily be deduced from the be- 
havior of the mean-field current shown in Fig. 4 of Ref. 1 1 
(see also Fig. 2 in Ref. 12). 

Including thermal as well as current-induced fluctuations 
by solving for the stationary solution of the Fokker-Planck 
equation (11) numerically and computing the shot noise with 
the help of Eq. (Al) leads to qualitatively the same effects as 
for the average current. This is exemplified in Fig. 12, where 
fluctuations lead to a softening of the borders in the v-v g plane 
delimiting the regions with a finite shot noise S s h — e 2 T/4. 



Appendix A: Shot noise 



Appendix B: Transition rates 



In this appendix, we briefly comment on how the shot noise 
properties of the setup considered in this paper (see Fig. 1) 
are influenced by the coupling (4) between the charge on the 
quantum dot and the mechanical degree of freedom. 

Within the adiabatic approach presented in Sec. II, the shot 
noise reads 

S Bh (ft)= f dxdpV st (x,p)S sh (n,x), (Al) 

where S s h(fl,x) is the quasistationary shot noise for fixed 
position x. Its zero-frequency limit (il -C T) 34 reads in the 
sequential-tunneling regime (hT -C k^T) 



S sh {0,x) 



4 



[/i» + /r(z)] [2 - h(x) - f R (x)} 

(A2) 

where /l/r(^) are the Fermi factors defined in Eq. (8). 

In the strictly adiabatic limit ujq/T —> 0, the current- 
induced fluctuation and damping coefficients are both vanish- 
ing [see Eqs. (9) and (10)], such that the stationary probabil- 
ity distribution function corresponding to the Fokker-Planck 
equation (11) is a Boltzmann distribution at temperature T. 
Assuming zero temperature, the (zero-frequency) shot noise 
(Al) reduces to S s h(0) = 6> s h(0,a; m ), with x m the global 
minimum of the effective potential (15). 11 Thus, the shot 
noise observes the same behavior as the mean-field current 
discussed in Ref. 11: In the v-v g plane, the shot noise has 
the same structure as the Coulomb diamonds delimited by 
slopes v ~ ±2%, the apex of these diamonds being given 
by Eqs. (16) and (17). The zero-frequency shot noise thus 



The transition rates Ti n and r out entering our approximate 
expressions for the average current (21) and the noise power 
spectrum (22) can be calculated using Kramers reaction rate 
theory. 32 ' 40 ' 41 In the weak damping regime, which is the exper- 
imentally relevant one for carbon nanotube-based resonators 
that can present very high quality factors, 9,10 the escape rate 
from the minimum located at X{ is given by 



7C 2nf 



i exp 



Vb,i 

f 



(Bl) 



where So(vb,i) is the abbreviated action at the barrier top, 
whose energy, seen from Xi, is denoted by Vh.i- Notice that 
Eq. (Bl) is valid as long as T <C «b,i an d in the weak damp- 
ing regime, i.e., j c So{vb,i) "C T. 41 With Eq. (Bl), and taking 
into account the probability that the system thermalizes in the 
well it jumps to, 41 we find for the transition rates 



r in = -y c u} —~h —— exp 



4T V2A„ 



Tout = 7cW — ~h [ | exp 



with 



2T V2A„ 



h(z) 




z 2 (l-z) 2 



z 2 + 2(l-z) 



(B2a) 



, (B2b) 



(B3) 
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